cpi uiaii 



O 

o 

C 

c3 



s 

a 

o 

o 



Pairing of Cooper pairs in a Josephson junction network contain- 
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PACS 74.20.Mn - Nonconventional mechanisms (spin fluctuations, etc.) 
PACS 74. 81. Fa - Josephson junction arrays and wire networks 
PACS 75 . 10 . Pq - Spin chain models 

Abstract. - We show how to induce pairing of Cooper pairs (and, thus, 4e superconductivity) 
as a result of local embedding of a quantum impurity in a Josephson network fabricable with 
conventional junctions. We find that a boundary double Sine-Gordon model provides an accurate 
description of the dc Josephson current patterns, as well as of the stable phases accessible to 
the network. We point out that tunneling of pairs of Cooper pairs is robust against quantum 
fluctuations, as a consequence of the time reversal invariance, arising when the central region of 
the network is pierced by a dimensionless magnetic flux ip = n. We find that, for <p> = 71", a stable 
attractive finite coupling fixed point emerges and point out its relevance for engineering a two 
level quantum system with enhanced coherence. 
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^-H , There is a large number of physical systems that can be 
■ mapped onto quantum impurity models in one dimension 
| PQ. Embedding a quantum impurity in a condensed mat- 
00 ter system may alter its responses to external perturba- 
' tions [2], and/or induce the emergence of non Fermi liquid, 
strongly correlated phases [3]. In quantum devices with 
T". ] tunable parameters impurities may be realized by means 
^ ■ of point contacts, of constrictions, or by the crossing of 
k>( [quantum wires or Josephson junction chains [JHZ]- For 
*j_J ■ instance, novel quantum behaviors have been recently ev- 
C$ ' idenced in the analysis of Y-] unctions of quantum wires [5] 
and of Josephson junction (J J) chains |7|. Here, we show 
how embedding a pertinent impurity in a JJ-chain may 
lead to the emergence of nontrivial symmetry protected 
quantum phases associated [8UTTj with the emergence of 
4e superconductivity in the network. 

While a standard perturbative approach works fine 
when impurities are weakly coupled to the other modes 
of the system (the "environment" ) , there are situations in 
which the impurities are strongly coupled to the environ- 
ment, affecting its behavior through a change of boundary 
conditions: when this happens, it is impossible to disen- 
tangle the impurity from the rest of the system, the per- 
turbative approach breaks down, and, consequently, one 



has to resort to nonperturbative methods, to study the 
system and the impurity as a whole. Such nonperturba- 
tive tools are naturally provided by boundary field theories 
(BFT) PQ: BFTs allow for deriving exact, nonperturbative 
informations from simple, prototypical models which, in 
many instances, provide an accurate description of exper- 
iments on realistic low dimensional systems ITS?] . In par- 
ticular, BFTs have been succesfully used to describe the 
dc Josephson current pattern in Josephson devices, such 
as chains with a weak link [Tgim] , and SQUIDs [151115] ■ 

In this letter, we analyze a Josephson junction network 
( J JN) , whose BFT description is given by a boundary dou- 
ble Sine-Gordon model (DSGM) [17] • The device is made 
by two half J J chains, joined through a weak link to a 
central Aharonov-Bohm cage C [18] , pierced by a (dimen- 
sionless) flux ip = $/$o($o = hc/(2e)) (see FigO}. For 
simplicity, we connect the outer ends of the chains to two 
bulk superconductors at fixed phase difference a. 

As we shall see, the dc Josephson current across the 
JJN is a periodic function of a of period 2tt, given 
by Ij{a) = E$(<p) sm(a) + 2E$ (ip) sm(2a). Varying 
ip changes the ratio E^(ip)/E^(ip); in particular, for 
tp = 7T, E^'^tt) = 0, and, thus, the period of Ij(ct) is 
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Fig. 1: The network. 

halved, signalling the tunneling of pairs of Cooper pairs 
(PCP)s through the JJN. A semiclassical analysis [H] 
already accounts for the two harmonics contributing to 
Ij{a): though, for a generic <p, single Cooper pair (SCP) 
tunneling is the dominating process for charge transport 
across the JJN (see top panel of FigJ^J) , for <p = tt, disrup- 
tive interference across C forbids SCP tunneling, and only 
allows for PCP tunneling, thus letting 4e superconductiv- 
ity to emerge in the JJN (see bottom panel of Figj2j). 

It is well known [15] that semiclassical and/or mean field 
approaches break down in one-dimensional JJNs, when the 
quantum phase fluctuations diverge logarithmically with 
the length of the system thus inducing a nonperturbative 
renormalization of the Josephson couplings. Here, we pro- 
vide a full quantum treatment of the JJN in Fig|TJ based 
on the BFT approach: We shall derive the dc Josephson 
current patterns and evidence the emergence- for a suit- 
able choice of the control and constructive parameters of 
the JJN- of a robust 4e superconductivity associated to a 
new attractive fixed point in the phase diagram accessible 
to the network. 

The central region is described by 
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and Xj,j+i = Xj - Xj+i + %■ 
Xj is the phase of the superconducting order parameter at 
the island j, V g is a gate voltage applied to each supercon- 
ducting island, J is the Josephson energy of each junction, 
E c is the charging energy at each site, and J z accounts for 
Coulomb repulsion between charges on nearest neighbor- 
ing junctions. In the charging regime, that is, -g- <C 1, 
and tuning V g so that V g = AZ+i, with integer A/", one may 
write He as a spin-1/2 XXZ-model [13], whose Hamilto- 
nian is given by He = - J Yl, 3 =o{ e% ^ Sj +l + h.c.} + 

J z J2j=o SjSj+i- Using the standard bosonization ap- 
proach [T2], the JJN Hamiltonian may be described by 
two spinless Luttinger liquids (LL), interacting with iso- 
lated spin-1/2 variables on C, with Hamiltonian given by 

ffjjN = H LL [{<P a }]+H T [{<I> a }, {8 a }; Si, S 3 ]+ffjr[Si, S 3 ], 

(2) 

with the Luttinger liquid Hamiltonian i/iLli'J'a}] = 
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Fig. 2: Leading tunneling processes across C: Top panel: Sin- 
gle Cooper pair tunneling for ip ^ n; Bottom panel: Tunnel- 
ing of pairs of Cooper pairs for tp — tt. 
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Hamiltonian H T [{$ a }, {6 a }; Si, S 3 ] 

+ h.c] 
Sf ( e L = l,e R = -1) 
+ St S3] + K z Si Sf. S 1; S 3 



and 



H K [S 1: S 3 ] - K[StS 3 + b{ s 3 j + K z ^a 3 . Si,» 3 are 
defined as in Figfl] <E> a is the collective plasmon field 
of the half chain a , while Q a (x,i) is its dual field, that 
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LL parameters are defined as g = tt/[2(tt - 



and £^££> = 12^11. The 
arccos( : f))], 

u = v / f(^/l-(f) 2 )/(arccos(f)), with A - E*/Ej, 
Vf — 2aEj, where Ej and E z are the Josephson energy 
and the Coulomb repulsion energy of the half chains, 
and a is the lattice step [20]. In deriving Eq.([2|), it is 
assumed that J/Ej < 1 and J z /E z < 1 (i.e., that 
C is weakly coupled to the chains); this allows to 
use Neumann boundary conditions at a: = 0, that is 

8$l(0) = d$fl(0) = q 
dx dx 

The couplings K,K Z in H T [$ L , $r, Q L , e R ; S u S 3 are 
dynamically generated by the interaction between C and 
the half chains, as it may be easily inferred from the renor- 
malization group (RG) equations for the dimensionless 
couplings JC = LK,JC Z = LK Z , and J = L 1 ~^J,J Z = 
J z . Indeed, employing standard BFT techniques [2"T1|2"2"] . 
one obtains „ dK , ; = K + [1 + cos(<p)] (J) 2 — : - — - 
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~ 0, show- 



ing that both JC and JC Z are dynamically generated when- 
ever ip =/= TT. 

Integrating the RG equations one sees that J(L) ~ 
L 1_ 2i, while JC(L) ~ L 2 ~s, that is, for ip 7^ tt, JC is al- 
ways more relevant than J . At variance, for ip = tt, no 
AT-coupling is generated and J is the only relevant cou- 
pling strength. For g > 1/2, the BFT description of the 
JJN allows to make very general statements regarding the 
regimes accessible to a network of finite size L. Namely, 
there will be a perturbative weak coupling regime, accessi- 
ble for small JC and J , and a non-perturbative strong cou- 
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pling regime, accessible when K, or J becomes ~ 1. Most 

importantly, there will be a renormalization group invari- 

_ i 

ant length scale ~ J 2 »- 1 , such that for L < L*, the 
JJN is in the perturbative weak coupling regime, while it is 
in the nonperturbative strong coupling regime for L > L*. 

To account for the last term of Eq. ([2]) , one may resort 
to a Schrieffer- Wolff transformation [23], which amounts 
to projecting over the ground state of C, by summing over 
its high energy states. At cp = tt, the ground state of C 
is twofold degenerate, while such a degeneracy disappears 
at cp ^ tt. As a result, for cp ^ tt, the central region C 
is effectively described by the boundary Hamiltonian Hb , 
given by 

H B [<p] = -E$(cp)co S [$-(0)]-E$(cp) C os[2<f>_(0)}, (3) 




2J Z 



2>K+K* + sin (f ) COS(^) K ( K+K zyi, 

while $_ = [$ L - $ R ]/V2. 

(Ej) 4 



with E$(cp) = cos(- 

At variance, for ip = n, H B [n] ~ cos[2$_(0)] 

I j (a) may be perturbatively computed as Ij(a) = 
lim^^oo —j 91 ^ , where Z is the partition function of the 

JJN, given by InZ sa \nZ Q -j^ cIt(Hb(t)), H b (t) being 
is the boundary interaction Hamiltonian in the (imaginary 
time) interaction representation, while Zq is the partition 
function at Hb = 0. From Eq.©, one obtains 
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Fig. 3: Josephson current vs. a for different values of ip and 
for J 2 /K 2 » 0.3: Top left panel cp = tt, Bottom left panel 
cp — 1.017T, Bottom right panel cp = 1.17T, Top right panel 

ip = 2n 



hi(uT/a), T being 



where the "instanton mass" M • 
the "instanton size", while V(x) 

Eyyif) cos(2a;). From the effective action in Eq.©, 
one gets the equation of motion MP 



TrE$(tp) sm[TrP] + 2%E$(ip) sin[27rP] = which, neglect- 
ing the "inductive term" cx 4, describes soliton solutions 
of the double-Sine Gordon model [17]. These are given 

with 



2L r 



Ij(a) = E$(<p) sm(a) + 2E$(cp) sin(2a) 



(4) 



with^'M = E$(cp) = byP(r) = Ea= ±1 |tan-i[ex P (^(ar + J R(^)) 



The ratio Eyy (cp)/Eyy (cp), measures the relative weight 
of SCP vs. PCP tunneling rate. One notices that, for 



R(cp) dehned by |sinh 2 



^ = tt, E$ (tt) = 0, while E$ (tt) f 0. Thus, at <p = tt 
PCP tunneling is the only allowed mechanism for charge 
transfer across C. In Fig.([3j), we plot Ij(a) vs. a for differ- 
ent values of cp. For <p ^ tt, Ij(a) is periodic, with period 
Aa = 2tt, while, for cp ~ tt, the period shrinks to Aa = tt. 

For g > 1, Hb is a relevant operator. Thus, when 
L/L* > f, the boundary Hamiltonian in Eq.© is 
the dominating potential term, and the field $_(0,r) 
takes values corresponding to minima of Hb', it obeys 
Dirichlet boundary conditions at both boundaries, yield- 
ing the mode expansion $_(a;, r) = a + (^p) ttP + 
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[14]. At the strongly inter- 



acting fixed point, instanton trajectories are the leading 
quantum fluctuations; they are described by imaginary 
time trajectories P(t), with P(t — > — oo) and P(r — > oo) 
corresponding to nearest neighboring minima of Hb- 
The instanton profile is derived from ^^pp = 0, where 

Ses[P] = JiUn^n] e" SE [*-] is obtained from the Eu- 
clidean action S£;[<i>_] for a spinless LL with parameters g 
and u, with a boundary interaction Hamiltonian given by 
Hb, after integrating over the oscilator modes of $-(£, r), 
{«„}. As a result, one gets 
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In Fig|U we plot Hb[$] vs. $ for various values of 
tp € [tt, 27r] and for -2tt < $(0) < 2tt. We see that, for 
cp =/= tt, the minima are separated by a distance 2tt. The 
corresponding instanton trajectories correspond to a "long 
jump" from P(— oo) = — 1 to P(oo) = 1 (represented by a 
solid arrow in Fig©. This may be regarded as a sequence 
of two "short instantons" (dashed arrows), separated by 
a distance 2R{cp) (see Fig©. As E^ (cp) becomes smaller 
(that is, as cp gets closer to tt), R(cp) increases, and even- 
tually diverges, as cp — > tt + . 

Though the degeneracy between the minima is broken 
by the 1/L-term, a pertinent tuning of the phase difference 
a allows to restore it. Indeed, for cp > tt, the degeneracy 
of the minima at $_(0) = 27rn±7r is restored by choosing 
a = 27m, n e Z , while, for cp < tt, one has to choose 
a = tt(1 + 2n) and, at cp = tt, a = 7r(i + 2n). 

In a BFT approach, long instantons are represented by 
the boundary vertex operators e ± 40 (° ! r ) j with scaling di- 
mension hi = g, while short instantons by the operators 
e ±§e(o : -r) ^ s Ca li n g dimension hi = |. For 1 < g < 4 
and cp = tt, short instantons are relevant perturbations of 
the strong coupling (Dirichlet) fixed point. As it happens 
for F-shaped JJNs, also here a new, time-reversal invari- 
ant, attractive finite coupling fixed point emerges in the 
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Fig. 4: Minima of the boundary potential Hs for differ- 
ent values of tp and for A*- 2 ) (ty) /A*- 4 ' (0) w 0.3 (in arbitrary 
units). The long (short) instanton trajectories are represented 
as solid(dashed) arrows: Top left panel tp = ir, Bottom left 
panel tp — 1.017T, Bottom right panel ip = 1.17T, Top right 
panel tp = 2n 

boundary phase diagram [7] , as a result of the twofold de- 
generacy of the ground state of C. At variance, for 1 < g 
and ip 7T, no finite coupling fixed point emerges, since, 
now, a possible departure from the Dirichlet fixed point 
should be due to long instantons, which are an irrelevant 
perturbation. For tp = n, the twofold degeneracy of the 
ground state of C is due to a Z2 symmetry of V[<fr_] man- 
ifesting its invariance under time reversal tp — > 2n — tp, 
while for tp 7^ n, time reversal is not anymore a symmetry 
of V[$_], as evidenced in Figs l4l5l Indeed, for tp = n + 6, 
two nearest neighboring minima at the top left panel of 
FigHJare splitted by an amount oc sin(|). As a result, for 
1 < g < 4 and for tp — n, the Z2 symmetry protects the 
robustness of PCP tunneling across C. 

We used quantum impurities as a resource for inducing 
local 4e superconductivity in a Josephson network, fab- 
ricable with conventional Josephson junctions. The pro- 
posed tunneling mechanism, realized for 1 < g < 4 and 
for tp = 7r, allows for the emergence of 4e superconductiv- 
ity, as a result of embedding an impurity in a Josephson 
network. Our analysis evidences that the emergence of 4e 
superconductivity is the signature of the presence in the 
dc Josephson current across the device of two distinct- and 
competing- periodicities, whose relative weight is tuned by 
the magnetic flux tp piercing the central region. Since, for 
a generic value of tp, the double sine-Gordon boundary po- 
tential has been normalized so that the Cooper pair charge 
2e is associated to the higher periodicity, one gets that, 
when - at tp equal to ir- the period is halved, the charge 
carriers should have charge 4e. Thus, the phenomenon 
analyzed in this paper is basically different from the mere 
phase difference renormalization taking place, for instance, 
in a series array of N equal Josephson junctions where, al- 
though the Josephson current is oc sin(a/iV), only a single 
harmonics is present. 
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Fig. 5: Profile of the instanton excitations -P(t) for different 
values of tp and for J 2 /K 2 « 0.3, M = 1: Top left panel 
tp — n, Bottom left panel tp = 1.017T, Bottom right panel 
tp = I.Itt, Top right panel tp — 2n 

Associated with 4e superconductivity there is, for 1 < 
g < 4, a new finite coupling attractive fixed point, which 
allows for the possibility of using the JJN as a two level 
quantum system with enhanced coherence [7J[24]. Fur- 
thermore, PCP tunneling is robust, as a consequence of 
the time reversal invariance, realized only for tp — n. The 
proposed mechanism exhibits intriguing similarities with 
the electron bunching phenomenon |25j . observed in shot 
noise measurements |26) on quantum dots in the Kondo 
regime. 
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